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We present magnetization (M) data of the d-metal alloy Nii_a;Va; at vanadium concentrations 
close to Xc ~ 11.4% where the onset of long-range ferromagnetic (FM) order is suppressed to 
zero temperature. Above Xc, the temperature (T) and magnetic field (H) dependencies of the 
magnetization are best described by simple nonuniversal power laws. The exponents of M / H ^ T~'' 
and M H" are related by 1 - 7 = a for wide temperature {lOK <T < 300K) and field {H < ST) 
ranges. 7 is strongly x dependent, decreasing from 1 at j: ~ a::^ to 7 < 0.1 for x=15%. This behavior 
is not compatible with either classical or quantum critical behavior in a clean 3D FM. Instead 
it closely follows the predictions for a quantum Griffiths phase associated with a quantum phase 
transition in a disordered metal. Deviations at the lowest temperatures hint at a freezing of large 
clusters and the onset of a cluster glass phase, presumably due to RKKY interactions in this alloy. 



PACS numbers: 71.27.+a, 75.40.-s, 75.50.Cc 

Magnetic phase transitions in metals continue to of- 
fer challenges to theory and experiment. In recent years, 
the focus has shifted from thermal transitions (such as 
the onset of ferromagnetism in nickel at a temperature 
of 630K [T|) to quantum phase transitions (QPTs) [2] 
that occur at zero temperature when a parameter such 
as pressure or chemical composition is varied. Spin fluc- 
tuations associated with continuous QPTs or quantum 
critical points (QCPs) are believed to be responsible for 
a variety of exotic phenomena including deviations from 
the fundamental Fermi liquid behavior of normal metals. 

Ferromagnetic and antiferromagnetic QCPs have been 
observed in transition metal alloys and heavy-fermion 
compounds (see Ref. ^ for a review). Quantum critical 
behavior is signified by singularities in thermodynamic 
and transport properties. According to the standard the- 
ory of ferromagnetic quantum criticality in 3D metals |4j , 
specific heat C, magnetic susceptibility x a-nd electrical 
resistivity p should behave as x ~ 2^-4/3^ C/T ^ log(T) 
and p{T) ^ j^s/s .^^j^gj^ approaching the QCP at low 
temperatures T. This was observed in Ni^Pdi_a; with 
X = 0.025, at least over a limited temperature regime 
[5]. However, most "clean" weak ferromagnets like MnSi 
[5], ZrZn2 [7] or NisAl [S] show deviations from the above 
predictions. The QPT becomes first order ^9j and is often 
accompanied by the appearance of novel phases. 

Many ferromagnetic binary alloys such as Nii_2:Cu2: 
or Nii_^Vj: in which T^, can be tuned by chemical sub- 
stitution X show a still more complicated behavior, even 
in the paramagnetic phase. In early investigations llOj. 
the existence of large magnetic clusters with giant local 
moments was proposed to describe the magnetization AI 
data of these inhomogeneous systems. 

Recent theories address the impact of disorder on 
QPTs more systematically (for a review, see Ref. 
We are interested in the case of metallic (itinerant) 
Heisenberg magnets. For these systems a strong-disorder 
renormalization group |12j predicts an exotic infinite- 
randomness QCP, accompanied by quantum Griffiths sin- 



gularities [13]. At such a QCP, thermodynamic observ- 
ables are expected to be singular not just at criticality 
but in a finite region around the QCP called the Grif- 
fiths phase. This is caused by rare spatial regions that 
are locally in the magnetic phase while the bulk is still 
nonmagnetic. The probability w of finding such regions 
is exponentially small in their volume V, w ^ exp(— 61^) 
with b a constant that depends on the disorder strength. 
Importantly, the characteristic energy scale e of a locally 
ordered region also depends exponentially on its volume, 
e ^ exp(— cy). Combining these two exponentials yields 
an energy spectrum P(e) ^ e'^^^. The nonuniversal Grif- 
fiths exponent X — b/c takes the value at the quantum 
critical point and increases with distance from criticality. 
This power-law spectrum gives rise to power-law quan- 
tum Griffiths singularities of many observables, including 
specific heat C ^ T^, susceptibility x ^ T^~^ and the 
zero temperature magnetization- field curve, M ^ . 

Quantum Griffiths singularities have attracted a lot 
of attention, but clearcut experimental verifications have 
been slow to arrive. Many heavy fermion compounds dis- 
play anomalous power-laws in C{T) and x{T^) |14j : and 
quantum Griffiths behavior was suggested as an expla- 
nation [TS]. However, in most of these systems, a sys- 
tematic variation of the exponents in accordance with 
theory could not be observed. Only recently, a partial 
confirmation could be found at the ferromagnetic QPT 
of CePdi_a;Rha; [IB]. It must also be noted that the in- 
terpretation of experiments in heavy-fermion compounds 
suffers from additional complications due to the Kondo 
effect which plays a crucial role for the magnetic prop- 
erties. It is thus highly desirable to observe quantum 
Griffiths singularities in a simpler system. 

In this Letter we therefore study the transition metal 
alloy Nii_^Va: as an example of an itinerant fcrromagnet 
in which Tc can be tuned to zero by chemical substitution 
while introducing strong disorder. We show that magne- 
tization and susceptibility close to the critical vanadium 
concentration indeed follow the quantum Griffiths see- 
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FIG. 1: Temperature-concentration phase diagram of 
Nii_2,Va; showing ferromagnetic (FM), paramagnetic (PM), 
quantum Griffiths (GP), and cluster glass (CG) phases. Cir- 
cles mark Tc found via the Arrott analysis. The grey line is an 
extrapolation of Tc{x). Also shown are Tmax defined by low- 
field maxima in x{T) and Tcross below which frozen clusters 
dominate x{T), leading to superparamagnetism. Inset: sat- 
uration magnetization ^sat vs x (determined as M(T < 5K, 
H > IT)). Data from Ref. [T^ for a; < 11% are included. 



nario over a wide temperature and magnetic field range. 

Polycrystalline spherical Nii_a;Va; samples with x = 
9 — 15% were prepared by arc melting from high purity 
elements (Ni 99.995%,V 99.8%) and annealed at 900 - 
1050° C. X-ray diffraction confirmed a single phase fee- 
structure with lattice constant a — (0.352 -I- 0.023x)nTO. 
Magnetization measurements were performed in a Quan- 
tum Design SQUID magnetometer from 1.8K - BOOK and 
magnetic fields up to 5T. The ac-susceptibility was mea- 
sured in a pick up coil in a dilution refrigerator down to 
0.05K and calibrated through the overlap with the mag- 
netometer data. All data shown are demagnetized. 

It is known that of Nii_a.Va; is rapidly reduced 
with increasing V-concentration x |T7j. As explained 
by Friedel [TB], Nii-x^x resides on a side branch of the 
Slater-Pauling curve: a V impurity (with 5 fewer elec- 
trons than Ni) creates a localized charge and a spin re- 
duction on the neighboring Ni-sites. This reduces the 
average spin moment by bfis/V from O.G/is/Ni leading 
to a critical concentration of about 12%. It also cre- 
ates large defects yielding an inhomogeneous magnetiza- 
tion density which makes Nii_a;Va; an ideal compound to 
study a QPT with significant "disorder" . 

The phase diagram resulting from our measurements 
is shown in Fig. [T] To get an overview, we first perform 
a standard Arrott analysis. Samples with x < 11% show 
mean- field behavior, i.e., parallel isotherms of the form 
H/M — a + bM^ as is common for itinerant magnets. 
(Mean-field behavior is expected outside the actual crit- 
ical region; if the standard theory [4] applied, it would 
describe the entire transition up to log. corrections.) Tc 
is extracted from the Arrott plots via the mean-field T- 
dependencies of magnetization and susceptibility. 
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FIG. 2: ac-susceptibility Xac of the x = 12.25% sample 
for several dc magnetic fields H and frequencies u of the (in 
phase) ac-magnetic field {Hac ~ O.IG) vs temperature T. In- 
set shows the frequency dependence of the maximum in x{T) 
Tmaa:- The fine follows dT^ax / d{\ogv) = 0.018K/dec(!/). 



For X > 11%, clear deviations occur from mean-field 
behavior (linear Arrott plots) [12] ; and the determination 
of Tc becomes sensitive to model assumptions. In addi- 
tion to Arrott plots, we analyze the differential suscepti- 
bihty x(T) = dM{T)/dH. It exhibits a field-dependent 
maximum at T^naxiH), indicating spin ordering or freez- 
ing. Figure [l] shows Tmax{H — >■ 0) estimated by a lin- 
ear extrapolation of the data taken at 0.5T to O.IT. 
Tmax{H — 0) is somewhat lower than Tc derived from 
the high field mean field analysis. We note that the mag- 
netization for X > 11% and higher fields H > 0.5T can be 
described by parallel isotherms in a modified Arrott plot 
(20| ((ff/M)(i/T) = a-h6M(i//9)) with exponents /3 = 0.5 
and an x-dependent j^x) < 1 |19| . However, this yields 
a finite Tc well above Tmax for all x < 15%. 

Within our error bars, Tmax{H — > 0) is definitely 
nonzero for x < 11% but zero for a; > 13%. To un- 
derstand the behavior at intermediate concentrations, 
we measure the ac-susceptibility Xac in ^ small ac-field 
Hac ~ O.IG down to 50mK as shown in Fig. [2] for the 
sample with x = 12.25%. A maximum in Xac{T) at 
Tmax = 0.19X marks spin freezing in OmT at a low fre- 
quency of = 340IIz. It is rapidly suppressed in small 
dc-fields and shifted to higher T . T„iax is dependent on 
V like in a spin glass, signifying irreversibility in this sys- 
tem. At higher T, no significant hysteresis in M{H) was 
found for all samples with a remanent field larger than 
the rest field of the magnet of the order of lOG. 

We emphasize that deviations from linear Arrott plots 
and the sensitivity of Tc towards the extrapolation pro- 
cedure already point to an unconventional QPT. More- 
over, the x-dependence of Tmax {H — >■ 0) in the accessible 
temperature region is better described by an exponential 
rather than a power law, making the determination of 
the Xc from finite-temperature data difficult. 

We now turn to the paramagnetic phase above the crit- 
ical concentration Xc ~ 11%. In the past, the susceptibil- 
ity at higher temperatures (T > 40 K) has been described 
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FIG. 3: (a) Low-field susceptibifity Xm = M/H — Xorb of 
Nii-^cVo; vs temperature T and (b) low-temperature magne- 
tization Mm ~ M—XorbH VS magnetic field H {x = 11 — 15%). 
Dotted lines indicate power laws for T > 10 K and H > 3000 
G in (a) and (b), respectively. Solid lines (shown for x > 12%) 
represent fits to the model ([2| with frozen clusters. 



[To] by a Curie- Weiss law x = C/(T - 6*) + Xorb, but this 
only works if the orbital contribution Xorb is allowed to 
vary by a factor of 3 with concentration x. In our anal- 
ysis, we keep the orbital susceptibility a;-independent at 
Xorb = 6 X 10~^eniu/niol which is the best fit parameter 
for X < Xc and lies within the reported estimates [21 . 

The resulting low- field {H = 100 G) spin susceptibility 
Xm ~ M/H — Xorb is showu in Fig. [sji for samples with 
X = 11 — 15%. At temperatures above 10 K, simple 
power laws describes the data well. We parametrize the 
power law by Xm{T) = A{kBT)~'^ . The exponent 7 is 
determined from fits between 30K and 300K (300K may 
seem a very high temperature for analyzing a QPT, but 
it is still well below the Curie point of Nickel at 630K; 
this high bare temperature scale is another advantage of 
our material.) Figure [4^ shows that 7 varies from about 
1 for X = 11.4% towards 0.03 for x = 15%. We also 
analyze the magnetization-field curve. Figure |3]d shows 
Mm = M — XorbH as a function of field at the lowest 
T — 2K for samples with x > 11%. For fields above 
H > 3000G « fesT/lO/is, M{H) follows a power law 
with an exponent a. Its value (shown in Fig. [4^) matches 
the susceptibility exponent, 1 — a = 7 for x > 12%. 

The results for Xm{T) and M„i{H) are in excellent 
agreement with the predictions for a quantum Griffiths 
phase outlined in the introduction if we identify the Grif- 
fiths exponent via A = a = 1 — 7. The critical concentra- 
tion Xc can now be determined from the condition A = 
which gives Xc ~ 11.4% using the susceptibility data. 
Right at criticality, the theory [T^] predicts extra loga- 
rithmic corrections (which are notoriously hard to verify) 
to the power laws. Our fits did not noticeably improve 
by including logarithmic terms. 

To combine the temperature and field dependencies (at 
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FIG. 4: (a) Exponents 7, 1 — a obtained from Xm ~ and 
Mm ~ H" vs V-concentration x. The plot also shows '^iot 
derived from Xm{T) below 10 K and "idyn extracted from the 
fit to (pi. The a;-dependence of 7 can be fitted to a power law, 
(1 — 7) ~ {x — XcY^ (dotted line), as expected according to 
[12] (with vip ~ 0.42) or to an exponential, — ln7 ~ {x — Xc) 
(solid line), (b) H/T scaling plot for x — 12.25% showing 
data at several const. T from 2K - 300K and H = 500G. The 
line is a fit to the phenomenological form discussed below ([T| . 

fixed x) of the magnetization we suggest the scaling form 

M = YifiH/ksT) (1) 

where Y is the scaling function and /i is a scaling mo- 
ment. The corresponding scaling plot for x = 12.25% 
is shown in Fig. |4]d. All M{H,T) data for temperatures 
above 14 K collapse, confirming H/T scaling. The scal- 
ing function Y is well approximated by the phenomeno- 
logical form Y{z) = A' /{I + ;j-2)t/2 where A' = A/ ^T* 
is a constant. We have produced similar scaling plots 
for the other concentrations. The resulting exponent 
7 matches that obtained by a direct fit of x{T) for all 
X between 11.4% and 15%. The scaling moment /i in- 
creases from \hb at X = 15% to 12^^ at x = 11.4% 
demonstrating the growth of the typical cluster size with 
X — > Xc |22j . An analogous scaling form was used to 
describe the H/T scaling in heavy fermions [55]. It also 
gives the correct exponent for the nonlinear susceptibility 
(X3 = d{M/H)/d{H^) ~ T"-3 for T > H). 

Having established the quantum Griffiths phase, we 
now turn to its limits. At T < lOK, deviations are ob- 
served from the "Griffiths" power laws. For instance, 7 
increases by about 50% (see Fig. [4| . This exponent does 
not match 1 — a. We thus believe that the behavior below 
lOK deviates from the quantum Griffiths scenario either 
due to the Vanadium distribution not being perfectly sta- 
tistical or because the rare regions arc not independent. 

It was shown ^25j that RKKY interactions between the 
rare regions lead to a dynamical freezing of the largest 
clusters at low T and to the formation of a cluster glass 
at even lower T. To explore this possibility, we model 
the zero-field susceptibility 

Xm{T) = Xl + Xdvn = A^/{kBT) + A/{kBT)^'-- (2) 

as the sum of a Curie term xi (describing the frozen 
clusters) and a Griffiths term with an exponent 'jdyn- 
This model describes the data in Fig. [3^ over the entire 
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temperature region above Tmax- A similar model can be 
formulated for the magnetization-field data in Fig. [3|d. 

We define a crossover temperature Tcross as the tem- 
perature where xi exceeds Xdyn 

(see Fig. [T|. It can be 
regarded as the boundary of the Grifhths phase. For 
X ~ 12.25%, Tcross ~ 7K which is much higher than 
the cluster glass temperature 0.19K found via the ac- 
susceptibility. This leaves room for a significant super- 
paramagnetic regime where independent frozen clusters 
dominate. We note that superparamagnetism can arise 
in systems with Ising spin symmetry even without rare 
region interactions |26| . However, our system does not 
show any indications of reduced spin symmetry. 

Further analysis of the Curie term requires insight into 
the structure of the rare region interactions. In a purely 
percolative scenario, one expects Ai ^ \x — Xc\~'^'' with 
the percolation exponent 7c — 1.8. Our data indeed show 
a divergence of Ai with x ^ Xc but with an exponent of 
about 2.6. Alternatively, one can successfully model the 
Curie term as a contribution from a number of frozen 
clusters of fixed moment which increases from 2^b for 
X = 15% to 15fj,B for X — 12.07% as Xc is approached. A 
more detailed discussion will be published elsewhere. 

We emphasize that the theory [T^ [O] for quantum 
Griffiths effects in metals was originally developed for 
antiferromagnets rather than ferromagnets. In fcrromag- 
nets, mode-coupling effects produce an additional long- 
range interaction [3] which renders the disorder pertur- 
batively irrelevant at the clean QCP [27^. However, this 



does not preclude the existence of Griffiths singularities 
because they are nonperturbative degrees of freedom. In 
fact, it was recently shown [55] that the physics of inde- 
pendent rare regions in a ferromagnet is the same as in 
an antiferromagnet. We thus believe that the quantum 
Griffiths scenario is applicable to our system, at least 
above the crossover temperature where interactions be- 
tween rare regions become important i29j . 

In summary, we have presented magnetization and sus- 
ceptibility measurements of the transition metal alloy 
Nii_a;Va; close to the critical concentration for the onset 
of ferromagnetism. While the finite-temperature phase 
transition in the concentrated Ni regime {x < 11%) is 
well described by mean-field behavior, the diluted regime 
with low or vanishing cannot be described in terms of 
conventional critical behavior. Instead, the data follow 
the predictions of a quantum Griffiths phase associated 
with an infinite-randomness QCP over a wide temper- 
ature and field region. Previous specific heat [23] and 
transport data [TU] support this scenario via anomalous 
power laws in a wide concentration range (even though 
they were not discussed in terms of a Griffiths phase). 
Deviations at lower T hint at individual freezing of large 
clusters before the system enters a cluster glass phase. 
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